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ABSTRACT 

Based on spectroscopic constraints derived from NLTE line formation, we explore the likely 
range of stellar parameters (T e s and logg) for the hyper-metal-poor (HMP) star HE 1327—2326. 
Combining the constraints from Balmer line profiles and the Cai/n ionization equilibrium, a 
subgiant stage of evolution is indicated. This result is further supported by spectrophotometric 
observations of the Balmer jump. If a higher T c ff value was used (as favoured by some photometric 
calibrations), the spectroscopic analysis would indicate a turnoff-point stage of evolution. 

Using a stellar-structure code that treats the effects of atomic diffusion throughout the star 
in detail, we evolve a low-mass model star to reach the HR-diagram position of HE 1327—2326 
after roughly 13 Gyr. While the surface abundances are modified significantly (by more than 
1 dex for the case of uninhibited diffusion) , such corrections can not resolve the discrepancy 
between the abundance inferred from the non-detection of the Li I resonance line at 6707 A and 
the WMAP-based primordial lithium abundance. As there are numerous processes that can 
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destroy lithium, any cosmological interpretation of a lower-than-expected lithium abundance 
at the lowest metallicities will have to await sample sizes of unevolved hyper-metal-poor stars 
that are one order of magnitude larger. The situation remains equally inconclusive concerning 
atomic-diffusion corrections. Here, attempts have to be made to better constrain internal mixing 
processes, both observationally and by means of sophisticated modelling. With constraints on 
additional mixing processes taken from a recent globular-cluster study the likeliest scenario is 
that HE 1327— 2326's surface abundances have undergone mild depletion (of order 0.2 dex). 

Subject headings: diffusion — stars: evolution — stars: fundamental parameters — stars: abundances 
— stars: Population II — stars: individual (HE 1327-2326) 



1. Introduction 

The most metal-poor stars provide us with 
chemical information about the early Universe and 
the onset of the cosmic chemical evolution. This 
is based on the assumption that the surface abun- 
dances of these old stars indeed reflect the compo- 
sition of their birth cloud, and have not changed 
during the long lifetimes. In principle, there are 
several internal and external processes which may 
have an impact on the surface composition of old 
stars, such as internal pollution of the outer layers 
by nuclear burning products from the stellar in- 
terior (in the case the star is sufficiently evolved) 
or accretion of interstellar gas during ~ 13 Gyr of 
orbiting though the Galactic potential (the latter 
pro cess was recently sho wn to be rather inefficient, 
see 



Frebel et al 



(j2008bh ) 

HE 1327—2326 is a low-mass star with extreme 
properties: its excessively low metallici t y in te rms 

™ (120051 )) is 



of iron ([Fe/H]« -5.6, iFrebel et al 



accompanied by large (3 -4 dex) over-abundances 
of the elements C, N and O relative to iron, a 
property shared by the o ther two stars k n own i n 
this metallic i ty re gime (jChristlieb et al.l (|2002l ). 



Norris et al 



( 2007k What these stars can tell 
us about the properties of the first generation 
of stars (metal-free Population III) is still de- 
bated. The abundance pattern has been inter- 
preted to reflect the nucl eosynthetic yie l ds of a 



single Pop HI super nova (llwamoto et al . (2005); 



Mevnet et al. ( 20061 )). but see Venn and Lambert! 
(|2008l ) for an alternative scenario potentially ex- 
plaining the abundance peculiarities of hyper- 
metal-poor (HMP) stars. 

Apart from these interpretative uncertainties, 
our ability to read the chemical fingerprints en- 
coded in the spectra of these stars is limited 
by the technical challenge to properly model the 
line-formation process taking place in their atmo- 



spheres. For the most part, the chemical abun- 
dances have been inferred using one-dimensional 
(ID) hydrostatic model atmospheres in local ther- 
modynamic equilibrium (LTE) with rough (esti- 
mated) corrections for departures from LTE (so- 
called NLTE effects). In the case of calcium, the 
only element observed in two ionization stages, 
ID LTE leads to a clear mismatch for the abun- 
dances derived from Cai and Can: the abun- 
dance from the Ca I resonance line at 4226 A gives 
an abundance 0.57 dex below the abundance de- 



rived from the Can K line at 3934 A (Aokietal 



2006) . This may reflect the shortcomings of classi- 



cal ID LTE analyses when applied to HMP stars. 
Explorations using three-dimensional (3D) hydro- 
dynamic model atmospheres have modified the 
abundance ratios (bringing down the large over- 
abundances by typically 0.7 dex for the CNO ele- 
ments) , but have been unable to resol ve the above- 
mentio ned discrepancy (reported by IFrebel et al 



(|2008ah to be 0.46 dex in 1D-LTE and 0.47 dex in 
3D-LTE). 

The LTE assumption is particularly question- 
able when analysing (hyper-)metal-poor stars, as 
their atmospheres are UV-transparent leading to 
strong photo-ionization. Furthermore, electronic 
collisions are less efficient due to the reduced num- 
ber of electron donors like Fe and Mg. Both these 
effects will drive the statistical equilibrium away 
from LTE. It is thus worthwhile to check whether 
NLTE can solve the obvious mismatch for the case 
of calcium. 

HE 1327—2326 is either an evolved main- 
sequence (MS) or early subgiant (SG) star at 
T eff «6200K, i.e., HE 1327-2326 is a star on 



the Spite plateau of lithium (jSpite and Spite 
Il982h formed b y stars with T eS > 580 K an d 
[Fe/H] < -1.3 (jCharbonnel and Primasl l2005fh 
Despite the uncertainty in the evolutionary status 
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of HE 1327—2326, the relatively unevolved nature 
of this old star makes it an ideal target for a Li 
measurement. Other unevolved metal-poor stars 
have been used to infer the primordial Li abun- 
dance. However, the Li doublet at 6707 A could 
not be detected in HE 1327-2326. The upper 
limit0 is log e(Li) < 0.7, in contrast to the Spite- 
plateau value of log e (Li) « 2.3 found in other 
stars. This finding currently lacks a satisfying 
explanation. 

In search for a solution, one can question the as- 
sumption that the surface abundances of old stars 
are constant over time. Even for single MS or SG 
stars, processes may be at work which systemati- 
cally alter the chemical surface composition. For 
solar-type stars, stellar-structure models allowing 
the atoms to move freely inside the star under the 
prevailing forces predict sizeable effects in particu- 
lar for old metal-poor stars. Such processes (grav- 
itational settling and radiative levitation, to name 
two prominent ones) are collectively referred to as 
atomic diffusion and require billions of years to al- 
ter the abundances of solar- type stars. Convection 
can suppress their efficiency, but atomic diffusion 
seems capable of depleting the thin low-mass con- 
vective e nvelopes of metal-po or stars appreciably 
(see e.g. Richard et al. (2002)), unless extra mix- 
ing processes stir up the boundary layer between 
the convective and the radiative zone. 

Observational evidence in favour of atomic 
diffusion in the presenc e of some extra mixing 



was recently pr esented (jKorn et al. 2006, l2007t 



iLind et al]l2008h for stars in NGC 6397 ([Fe/H 
—2). Based on observed abundance trends be- 
tween stars from the turn-off point (TOP) to 
the red-giant branch (RGB), the mixing efficiency 
needed in the stellar-structure models was empir- 
ically determined. This mixing falls in the range 
expected from obs ervational constrain ts set by 
the Spite plateau ([Richard et al.l 120051 ): surface- 



abundance alterations due to atomic diffusion are 
then limited to ±0.3dex and help to harmonize 
the inferred lithium abundances with the predic- 
tions from WM AP-based Big-Ban g nucleosynthe- 



sis calculations ( Korn et al. 2006h 



In this article, we thus explore to which extent 
the surface abundances of HE 1327—2326 may be 



affected by atomic diffusion. After describing the 
observational material (Sect. 2) we derive stellar 
parameters using spectrophotometry and NLTE 
diagnostics (Sect. 3). Based on these, we present 
upper limits on the abundance variations induced 
by atomic diffusion and discuss their consequences 
(Sect. 4). In Sect. 5, we summarize our conclu- 
sions. 

2. Observations 

We combine data obtained in 2004 and 2005 
with Subaru/HDS and VLT/UVES, respectively. 
The U VES spectra (R = 40 000, iFrebel et al ' 



( 20061 )) are used for the Balmer lines (Ha and 
H/3) and the Ca IR triplet line s . The HDS spec- 



tra (R = 60 000, lAoki et alj (|2006h ) cover Ca 



K 3933 A and the Cai resonance line 4226 A. In 
all cases, the signal-to- noise ratio exceeds 200:1. 
The limiting factor for an abundance analysis of 
HE 1327—2326 is thus not to be sought in the 
quality of the input data. 

Spectrophotometric observations of HE 1327—2326 
were obtained with the Double Beam Spectro- 
graph (DBS) on the 2.3m telescope at Siding 
Springs Observatory, Australia. In June 2006, 
a 600B' spectrum from 3200 A to 5200 A was ob- 
tained. In January 2009, a 300B' spectrum was 
taken covering the wavelength range from 3200 A 
to 6200 A. The spectra were taken with the atmo- 
spheric dispersion along the slit to ensure that the 
relative absolute fluxes were accurately measured. 
The spectra were corrected for continuous atmo- 
spheric extinction and wavelength calibrated using 
Ne-Ar arc exposures. They were also divided by 
a continuous white dwarf spectrum to remove the 
sp ectrogra p h and dichr oic response as described 
by iBesselll (|l999[ 120071) . The spectra were then 
scrunched and flux corrected using more than five 
spectrophotometric standards selected each night 
from the Next Generation Stellar Library STIS 
observations of Michael Gregg as re- reduced by 
Heap and Lindler . The relative absolute fluxes 
of the NGSL stars are precise to within two per- 
cent. The 300B and 600B spectra show excellent 
agreement. 



lr The customary logarithmic 12-scale is used here: 
log e(X) = log n x /n H + 12 



2 http: / /www. stsci.edu/hst/HST_overview/documents/ 

calworkshop / workshop2005 /papers/ gregg.pdf 
3 http: / /archive. stsci.edu/prepds/stisngsl/index. html 
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3. Stellar parameters 

A crucial ingredient in this investigation is the 
knowledge of HE 1327— 2326's evolutionary sta- 
tus. So far, it has not been possible to distin- 
guish whether HE 1327-2326 is an MS or an SG 
star. Luckily, this has no significant impact on 
the abundances derived from minority species (like 
Cai and Fei) and only a small one for majority 
ones (like Can, « 0.2 dex). For C, N and O, how- 
ever, there are differences of up to w 0.5 dex be- 
tween the two cases, as these abundances are ex- 
clusively derived from molecular features. Since 
the CNO abundances with respect to iron are ex- 
tremely high in this star, these differences have not 
greatly hampered the interpretation of the overall 
abundance pattern. But with increasing sophisti- 
cation of the modelling it is crucial to determine 
the evolutionary status of this star. This is also the 
case when comparing to predictions from stellar- 
structure models. 

3.1. Effective temperature 

From four broad-band indices calibrated on 



the in frared-flux method (IRFM) of lAlonso et al. 
(jl996h . a mean effec tive temperature of T P ff = 618 K 
was derived (jFrebel et al.1 120051 : lAoki et al.ll2006l ). 
Balmer lines in LTE seemed to favour a some- 
what lower effective temperature, T o ff«6000K, 
irrespectiv e of the broadenin g theory applied (see 
Table 7 of lAoki et al.1 (|2006h ). 



We have performed NLTE calculations for Hi 



using the method described in iMashonkina et al 
(|2008l ). In brief, the model atom includes lev- 



els with principal quantum numbers up to n = 19. 
Collisional rates include inelastic collisions with 
electrons and hydrogen atoms. For the latter, 
we us e the formula of ISteenbock and Holwegeil 
([19841 ) . Since it provides only an order of magni- 
tude estimate, we apply it with a scaling factor Sh 
= 2 as favoured by the consistency of the effective 
temperatures derived from two Balmer lines, Ha 
and H/3, in four metal-p oor ([Fe/H] < —2) stars 



(|Mashonkina et aLlExjj ). 



For the model atmospheres investigated in 
this study, the ground state and the n = 2 
level keep their thermodynamic-equilibrium level 
populations throughout the atmosphere. Fig. 
[T] shows the behaviour of the departure coeffi- 
/n^ of the relevant levels 



of H I as a function of continuum optical depth 
75000 referring to A = 5000A in the model with 
T cff /log 5 /[Fe/H] of 6180/3.7/-4.34. Here, nf LTE 



and 



.LTE 



are the statistical-equilibrium and ther- 



mal (Saha-Boltzmann) number densities, respec- 
tively. Departures from LTE for the n — 3 level 
are controlled by Ha. In the layers where the 
continuum optical depth drops below unity Ha 
serves as the pumping transition resulting in an 
over-population of the upper level. For instance, 
b(n — 3) ~ 1.08 around logr 50 oo = —0.5 where the 
core-to-wing transition is formed. Starting from 
logTsooo ~ — 3 and further out, photon losses from 
the wings of H Q cause an under-population of the 
n = 3 level. We obtain weak NLTE effects for the 
H/3 profile beyond the core (which is not consid- 
ered in the analysis), independent of the applied 
scaling factor for hydrogen collisions. For Ha, 
NLTE leads to a weakening of the core-to-wing 
transition compared to the LTE case. Exactly 
this part of profile is sensitive to T e g variations. 

Table 1 summarizes our T e g determination. In 
NLTE, both Ha and H/3 indicate an effective 
temperature close to T e ff = 6120K (see Figs. 2 
and 3) putting the new spectroscopic effective 
temperatures in much better agreement with the 
photometric estimates. For reasons of consis- 
tency with previous publications in this series, we 
adopt T e ff = 6180K as our current best estimate 
(a significantly higher effective temperature seems 
possible using different calibr ations of the IRFM 
(Ramirez and Melendez|[2005h and we briefly dis- 
cuss the impact of this scenario below). 

3.2. Surface gravity 

The surface gravity of hyper-metal-poor stars 
is, generally speaking, difficult to constrain ob- 
servationally. The Balmer jump can serve as 
a p hotometric su rf ace-gr avity indicator (see, 
e.g., lOnehag et ail (j2008l) V Figure 4 shows 
comparisons between the two DBS observa- 
tions and model fluxes of two MARCS mod- 



cients 6, 



.NLTE 



els (jGustafsson et al.1 120081 ) with T off = 6180K, 
[Fe/H] = —5 and logg between 3.5 and 3.75 
{E{B - V) = 0.075 was assumed). As the re- 
cent '300B' spectrum was obtained in bright time, 
higher weight is given to the '600B' observations. 
A spectrophotometric surface gravity of close to 
log g = 3.65 is derived. 

Ionization equilibria are a standard spectro- 
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Fig. 1. — Departure coefficients logb^ for the five lowest levels of Hi in the model atmosphere 6f 80/3.7/— 4.34. 
Tick marks indicate the locations of core-to-wing transition formation depths for H a and (0.85 to 0.95 
in normalized flux). 



Table 1: Spectroscopic effective temperatures as derived from Balmer lines in NLTE. 



star 


HckNLTE 


H/?NLTE 


error a 


HE 1327-2326 


6120K 


6120K 


±150 K 


Sun 


5780 K 


5780 K 


+20K, -60K 



"estimated total error budget including observational uncertainties 
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Fig. 2.— Fit to the observed Ha profile. An effective temperature of T e g = 6120K is indicated in NLTE. 
The predicted NLTE profile for T e g = 6180 K is indicated by the upper grey line, the LTE one by the lower 
one. 
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Fig. 3. — Fit to the observed H/3 profile (black line). An effective temperature of T c g — 6120 K is indicated 
in NLTE. The predicted NLTE profile for T cff = 6f80K is indicated by the grey line. 
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3500 4000 4500 5000 5500 6000 

Wavelength Angstroms 

Fig. 4. — Comparison between spectrophotometric observations of HE 1327—2326 and model-atmosphere 
predictions (T off = 6180 K, [Fe/H] = -5, E(B - V) = 0.075). The strength of the Balmer jump indicates a 
surface gravity of log g = 3.65. 
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Fig. 5. — Departure coefficients logbi for selected levels of Cai and Can in the model atmosphere 
6180/3.7/— 4.34. Tick marks indicate the locations of line-centre optical depth unity for Cai 4226 A, Can 
3933 A and Can 8498 A. 



9 



scopic method of constraining logg. The typi- 
cal attainable precision is 0.3 dex when line-to- 
line scatter is of order 0.1 dex. The accuracy 
can be increased by considering NLTE effects. 
The present investigation of Cai/ll in NLTE 



is based on our earlier work (Mashonkina ct al 



20071 ) where atomic data and the problems of 
calcium line formation were considered in de- 
tail. The model atom used is up-to-date and 
includes effective collision strengths for Ca n tran- 
sitions from r e cent R— matrix calculations of 
Melendez et al.l (|2007j ). Inelastic collisions with 
neutral H particles are accounted f o r usin g the 
formula of Ste enbock and Holwegeri (Il984i ) with 
a scaling factor Sh =0.1 as empirically deter- 
mined from their different influence on Cal and 
Call lines in Solar and stellar spectra. As in the 
case of H, such a scaling is warranted, as the 
Steenbock and Holweger ( 19841 ) recipe yields col- 
lisional rates with large systematic uncertainties. 
Non-LT E calculations are perf ormed for two 



MAFAGS (jFuhrmann et all 119971 ) model atmo- 



spheres with effective temperature T G ff = 6180K 
and iron abundance [Fe/H] = —4.34, but with dif- 
ferent surface gravities, log g = 3.7 and 4.5 (the in- 
put calcium abundance is iteratively adjusted to 
match the derived value). The behaviour of the 
departure coefficients and the line source functions 
is similar for both models. In the first place, we 
are interested in the behaviour of the lower and up- 
per levels from which the investigated lines, Cal 
4226 A, Can 3933 A and Can 8498 A, arise. Fig- 
ure 5 shows the departure coefficients of the se- 
lected levels as a function of continuum optical 
depth t 50 oo in the model 6180/3. 7/-4.34. 

For the minority species Cal, the main NLTE 
effect is over-ionization. It is caused by superther- 
mal radiation of non-local origin below the thresh- 
olds of the Ca I ground state 4s 2 X S (Athr = 2028 A) 
and the low-excitation levels 4p 3 P° ( A t hr = 2930 A) 
and 3d 3 D (Athr = 3450 A). The Ca I resonance 
line is weaker relative to its LTE strength not 
only due to the general over-ionization but also 
due to the larger departure coefficient of the 
upper level compared to that for the lower 
one resulting in a line source function Si u ~ 
fe(4p 1 P°)/6(4s 2 1 S)B U > B v and diminished line 
absorption. The NLTE abundance correction 
A NLTE = loge(NLTE) - log e (LTE) = +0.31 dex 
and +0.25 dex for the models 6180/3.7/-4.34 and 



6180/4.5/-4.34, respectively. 

In the line formation layers, Can dominates 
the element number density. No process seems 
to affect the Call ground state population and 
4s keeps its thermodynamic-equilibrium value. 
Call resonance- line pumping produces a slightly 
enhanced excitation of the level 4p in the at- 
mospheric layers between log T5000 = — 0.3 and 
— 1.3. In the very metal-poor atmosphere of 
HE 1327-2326, the metastable level 3d is only 
weakly coupled to the ground state and follows 4p 
until photon losses in the Can 8498 line start to 
become important. Departures from LTE for Can 
8498 are weak resulting in Anlte = — 0.02 dex 
and -0.01 dex for the models 6180/3.7/-4.34 
and 6180/4.5/— 4.34, respectively. Because of the 
larger line strength (and thus shallower forma- 
tion depth) of Can 3933, this line suffers from 
somewhat larger NLTE effects: —0.15 dex and 
—0.06 dex for the subgiant and dwarf model, re- 
spectively. 

We note that NLTE effects for Cai/n in 3D 
line formation calculations are currently unex- 
plored. For the minority species Li I in metal- 



(2003) 



poor stellar atmospheres, iBarklem et al. 
concluded that 3D + NLTE is not greatly different 
from ID + NLTE with ID + NLTE results overes- 
timating the Li abundance by up to 0.1 dex. Sim- 
ilar effects can be expected for Cai. For Can, the 
departures from LTE are likely to be small, similar 
to the ID case. 

Table 2 summarizes the results concerning cal- 
cium. Good agreement (within 0.05 dex) is ob- 
tained for the three different lines when NLTE and 
a surface gravity of log g = 3.7 are used (numeri- 
cally, a slightly lower log g value is indicated which 
we disregard for reasons explained in Sect. |4j. At 
the dwarf- like surface gravity, the log ^-sensitive 
line Ca 11 8498 is more than 0.4 dex discrepant from 
Cai 4226. From this and the line-to-line scatter 
we derive an uncertainty for log g of 0.2 dex. 

A rigourous NLTE modelling of Cai/n in ID 
is thus capable of fully accounting for the Cai/n 
abundance discrepancy seen in previous analyses. 
HE 1327—2326 is therefore a subgiant, a con- 
clusion in accord wi t h ind ications from Balmer 
lines fsee lAoki et~ai1 (|2006l )). The [Ca/Fe] ratio 
is +0.47 dex which is not an unusual value for a 
halo star. HE 1327—2326 is a benchmark case for 
future NLTE studies for Cai/n in 3D model at- 
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Table 2: Calcium abundances as derived from individual lines. Like in Aoki et al. ( 20061 ). a microturbulence 
value of 1.7km/s and 1.5km/s was assumed in the line formation calculations for the subgiant and dwarf, 
respectively. 



line 




loge(Ca) in NLTE 


loge(Ca) in LTE 


log E (Ca) in NLTE 


loge(Ca) in LTE 






@ log g = 3.7 


@ log g = 3.7 


@ log g = 4.5 


@ log g = 4.5 


3933 A 


128.9 mA 


1.17 


1.32 


1.29 


1.35 b 


4226 A 


2.7mA 


1.22 


0.91 


1.15 


0.90 


8498 A 


9.6 mA 


1.27 


1.29 


1.56 


1.57 


all lines ± 1(7 




1.17° ± 0.05 


1.17 ± 0.23 


1.33 d ± 0.21 


1.27 ± 0.34 



a as determined from the the theoretical spectra after profile fitting 

b In this analysis, the fit is based on the wings, as the wings and core cannot be fitted simultaneously. 
c corresponds to [Ca/Fe] =0.47 on the scale of lAoki et al. (2006) 
Corresponds to [Ca/Fe] =0.62 on the scale of lAoki et aD (2006) 



mospheres. 

We estimate the effects of crucial atomic data 
on the accuracy of the NLTE Ca abundances 
derived from Cai 4226 and Can 8498 (the two 
most critical lines for our purpose) using the 
model 6180/3.7/-4.34. The largest effects are 
found upon disregarding hydrogenic collisions al- 
together. But even in that case, none of the line 
abundances changes by more than 0.1 dex and our 
conclusion about HE 1327—2326 being a subgiant 
holds. 

4. Stellar-structure model predictions 
4.1. Model description 

The stellar evolutionary models used here are 
calculated with the physics of particle transport 
that is known from first principles. In radiative 
zones this implies taking into account atomic dif- 
fusion including gravitational settling, thermal dif- 
fusion and radiative accelerations, in addition to 
the purely diffusive term. The diffusion velocity of 
each species is computed with the equations as de- 
veloped in Burgers (1969) which take into account 
the interactions between all diffusing species. 

The detaile d treatment of atom ic diffusion 
is described in iTurcotte et alj (|1998l ). In these 
models, the Rosseland opacity and radiative ac- 
celerations are computed in each layer for the 
exact local chemical composition using OPAL 
monochromatic opacities and for each time step 
during the evol ution. The radiat ive acceler- 

(1 19981) with co r- 



Richer et al 



ations are from 
rections for redistribution from IGonzalez et al 



(1995) and lLeBlanc et all (2000). Convection and 
semi-convection are modelled as diff usion pro- 

(l2000l) 



cesses as described in Richer et al 



and 



Richard et alj (|200ll) . 

In one set of models, additional mixing is in- 
cluded in the diffusion equatio n by adding an effec- 



tive diffusion coefficient as in Richer et al. I (|2000l) 



and iRichard et al" ( 2001 ). The parameters spec- 
ifying this effective diffusion coefficients are indi- 
cated in the name assigned to t he model, wh i ch are 

(2007) 



the same as those used in the Korn et al 



study on NGC 6397. This additional mixing is in- 
troduced in an ad-hoc manner: here no attempt 
is made to connect it t o physical processe s like 
rotation. As shown by iRicher et al. ( 2000l ). the 
most important factor is the radial extent of the 
mixed layers. All models considered here were as- 
sumed to be chemically homogeneous on the pre- 
main sequence with Solar relative mass fractions 
as given in Table [3j To mimic the abundance pe- 
culiarities of HE 1327—2326, the relative concen- 
trations of the CNO elements were increased by 
+4 dex compared with the scaled Solar mixture. 
As the most abundant metals in HE 1327—2326, 
they have structural effects on the models. 

4.2. Best-fitting models 

We computed a variety of models with masses 
around 0.78 M© for this study, both purely dif- 
fusive ones and ones including additional mix- 
ing. The mass and initial composition of each 
model was varied to roughly meet the observa- 
tional constraints of HE 1327—2326: the surface 
abundance of iron should be within a factor of two 
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Fig. 6. — Position of HE 1327—2326 relative to evolutionary models constructed from the two diffusion 
models. 
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Table 3 

Initial chemical composition of the two different diffusion models. 
The abundance variations after « 13.3 Gyr can be estimated from 

Figure [7] 



Element 


Note 


Mass fraction 


[Fc/HJ = -5.61, T6.0 


[Fe/HJ = -7.11, diffusive 


H 




7.642X1CT 1 


7.645xl0~ 1 


4 He ... 




2.354xl0~ 1 


2.354X10" 1 


7 Li .... 




2.0xl0~ 9 


2.0X10 -9 


12 C ... 




8.633X10 -5 


2.763xl0~ 6 


N 




2.647xl0~ 5 


8.471xl0~ 7 


O 




2.409xl0~ 4 


7.708X10 -6 


Nc .... 




4.927xl0~ 9 


1.577xl0~ 10 


Na .... 




9.989X10" 11 


3.197xl0~ 12 


Mg .... 




1.875xl0~ 9 


6.002xl0~ 11 


Al 




1.623xl0~ 10 


5.194xl0~ 12 


Si 




2.023x10"° 


6.473xl0~ 41 


P 




1.748X10 -11 


5.594xl0~ 13 


S 




1.056X10 -9 


3.380X10 -11 


CI 




2.248X10 -11 


7.192X10 -13 


Ar .... 




2.697X10" 10 


8.631X10 -12 


K 




9.989X10 -12 


3.197xl0~ 13 


Ca .... 




1.873X10 -10 


5.994xl0~ 12 


Ti 




9.989X10 -12 


3.197xl0~ 13 


Cr .... 




4.994X10" 11 


1.598X10 -12 


Mn .... 




2.747xl0~ 11 


8.791xl0~ 13 


Fo 




3.586xl0~ 9 


1.147xl0~ 10 


Ni 




2.223xl0~ 10 


7.112xl0~ 12 


Z 




3.545xl0~ 4 


1.135xl0~ 5 



increased by +4 dex relative to the assumed solar mixture 
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Fig. 7. — Surface-abundance variations as predicted by the two diffusion models, relative to the initial 
abundances given in Table 3. 
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of [Fe/H] = -5.8 (the aver age of the ID and 3D 
LTE abundances de rived bv lAoki et al. ( 20061 ) and 
Frebel et all (|2008ah ) when the model star reaches 
the error ellipse of HE 1327—2326 on the subgiant 
branch after roughly 13Gyr (see Figure 6). 

Two models of 0.785 M fulfilled these a- 
posteriori requirements. The first one has an 
initial [Fe/H] of [Fe/H].; = — 7.11 and was com- 
puted with atomic diffusion and radiative acceler- 
ations, but without additional trans port processes 



i n the radiative zone. As shown by [Richard et al 



(2002) for higher metallicities (between —2.31 and 
—4.31), such a star is expected to display strong 
surface abundance variations close to and past 
the TOP. The second one ([Fe/H] = -5.61) addi- 
tionally includes additional transport as needed 
to explain abundance variations of NGC 6397 



(jKorn et al.ll2006L 120071 ) 



4.3. Birth-cloud vs. present-day abun- 
dances 

As can be seen from Figure 7, abundance cor- 
rections are large for the model with uninhibited 
diffusion: after w 13.4 Gyr the M = 0.785 M 
[Fe/H], = — 7.11 model has reached the subgiant 
branch and its surface lithium is depleted by al- 
most 1.2 dex. Simultaneously, radiative levita- 
tion has increased its convection-zone abundances 
of iron by 1.2 dex. Other elements, e.g. alu- 
minium and nickel, are affected by an even larger 
amount. We compared the diffusion-corrected 
abundance pattern of HE 1327—2326 to that of 
HE 0107-5240 and find that the abundance pat- 
terns of the two stars do not become more simi- 
lar. When additional mixing is introduced into the 
model, the surface-abundance changes are largely 
moderated: helium and lithium only settle by 
~ 0.3 dex (incidentall y, the same amou nt as pre- 



(|2007h ). iron by 



dieted for NGC 6397. lKorn et all 
only « 0.1 dex. 

If HE 1327—2326 turned out to be significantly 
hotter (T ff = 6450 K), then a surface gravity of 
log <?~ 4.0 would be derived in the framework of 
this analysis. From the stellar-structure point of 
view, this is a meaningful combination of stellar 
parameters (see Figure 6), albeit not supported 
by our understanding of the formation of Balmer 
lines. The same diffusion models could be used to 
predict corrections to the surface abundances and 
these would turn out to be quite similar within a 



given set of model assumptions. This is owed to 
the fact that the change in stellar parameters effec- 
tively moves the star along its fast early post-TOP 
evolutionary track. 

Given our present ignorance of the extent of the 
additional mixing (and how it scales with, e.g., 
metallicity), it is difficult to draw firm conclusions 
about which of the two diffusion scenarios is more 
realistic. It is conceivable that the mixing effi- 
ciency declines towards lower metallicities, as the 
convective envelope decreases in extent and mass. 
Nonetheless, even in the absence of mixing the cor- 
rection is not large enough to reconcile the up- 
per limit on the lithium abundance (log e (Li) < 
0.7) with the WMAP-based primordial lithium- 
abundance predictio n of Big-Bang nucle osynthe- 
sis (loge(Li) = 2.63. ISoergel et all (120071 ) ). If the 
low lithium content of HE 1327—2326 turned out 
to be the rule among unevolved HMP star, then 
this may tell us something specific about the en- 
vironment these stars were able to form in. More 
work is needed to identify similar stars in the outer 
halo. 

Until better constraints can be set on the mix- 
ing efficiency, we recommend to use the abundance 
corrections based on the model with extra mixing. 
These corrections are moderate and generally neg- 
ative. The birth cloud from which HE 1327-2326 
formed was somewhat more metal-rich than the 
surface abundances of HE 1327—2326 show at the 
present time. 

As discussed in iKorn et all (120071) ). helium dif- 
fusion has a structural effect on the atmosphere 
of a star and can be mapped as a shift in the 
surface-gravity scale. The gravitational settling 
of helium entails a shift of the surface gravity 
of « +0.05 dex. This correction was consid- 
ered when deriving our best current estimate of 
log ff = 3.7 ±0.2. 

5. Conclusions 

We have derived spectroscopic stellar parame- 
ters for HE 1327-2326 using ID NLTE techniques: 
Balmer lines profiles in NLTE yield T cff = 6120 
± 150 K confirming the photometrically deter- 
mined effective temperature of T e g = 6180K. 
The calcium ionization equilibrium in NLTE 
points towards a subgiant surface gravity of 
log<? = 3.7 ± 0.2, as does spectrophotometry of 
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the Balmer jump. Using an effective temperature 
of 6450 K (as favoured b y a re vision of the IRFM, 
Ramirez and Melendea (|2005[ )). results in a sur- 
face gravity typical for turn-off stars. While this 
is also a sensible stellar-parameter combination 
from the point of view of stellar evolution and the 
Cal/ll ionization equilibrium, it finds no support 
from Balmer-line analyses. 

The exploratory investigations into how much 
HE 1327— 2326's surface abundances may be af- 
fected by atomic diffusion show that large (more 
than 1 dex) corrections are in principle possi- 
ble. The material from which HE 1327-2326 
formed would then have been more metal-poor 
([Fe/H], w —7) and significantly more lithium- 
rich (loge(Li)i < 1.9). However, there is no 
compelling evidence in favour of uninhibited dif- 
fusion in stars that have convective envelopes. If 
additional mixing below the convective envelope 
is postulated (as n eeded to explain abundance 
trends in NGC 6397. lKorn et al.l (|2007l) ). then the 
abundance corrections are limited to a factor of 
two. Assuming the additional-mixing efficiency to 
be independent of metallicity, the corrections are 
thus moderate given the uncertainties associated 
with present-day 1D-(N)LTE or 3D-LTE abun- 
dances of HE 1327-2326 and other HMP stars. 

It is too early to draw firm conclusions about 
our ability to reliably read the chemical finger- 
prints in the atmospheres of stars like HE 1327—2326 
and their nucleosynthetic significance as such. 
More such stars need to be identified in deep wide- 
angle surveys to leave the realm of low-number 
statistics. Efforts in this direction are under way. 
Systematic comparisons between chemical abun- 
dances in HMP dwarfs, subgiants and giants may 
then reveal the efficiency of mixing moderating 
atomic diffusion. Simultaneously the modelling 
has to be improved. This will mean replacing the 
hydrostatic approximation with full radiation hy- 
drodynamics, both in stellar-structure and stellar- 
atmosphere modelling. A similar transition is re- 
quired for photospheric line-formation calculations 
in terms of statistical equilibria rather than LTE, 
as documented by the successes of this NLTE 
study. 
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